Renyi Entropy and Parity Oscillations of the Anisotropic Spin-s Heisenberg Chains in 

a Magnetic Field 
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Using the density matrix renormalization group, we investigate the Renyi entropy of the 
anisotropic spin-s Heisenberg chains in a z-magnetic field. We considered the half-odd-integer spin-s 
chains, with s = 1/2, 3/2 and 5/2, and periodic and open boundary conditions. In the case of the 
spin-1/2 chain we were able to obtain accurate estimates of the new parity exponents and 
that gives the power-law decay of the oscillations of the a— Renyi entropy for periodic and open 
boundary conditions, respectively. We confirm the relations of these exponents with the Luttinger 
parameter K, as proposed by Calabrese et al. [Phys. Rev. Lett. 104, 095701 (2010)]. Moreover, 
the predicted periodicity of the oscillating term was also observed for some non-zero values of the 
magnetization m. We show that for s > 1/2 the amplitudes of the oscillations are quite small. 



and get accurate estimates of p^' and p'^' become a challenge. Although our estimates of the new 
universal exponents and pL"' for the spin-3/2 chain are not so accurate, they are consistent with 
the theoretical predictions. 

PACS numbers: 03.67.Mn, 05.70Jk, 75.10.Jm 



I. INTRODUCTION 

The observation that entanglement may play an im- 
portant role at a quantum phase transition has moti- 
vated many studies on the characterization of the critical 
phenomena by using quantum information concepts ji^— 
Quantum spin chains have been proven as useful labora- 
tories to investigate the interconnection of entanglement 
and quantum criticalityi^"— Although does not exist yet 
an universal measure that quantifies the entanglement;^ 
the von Neumann entropy and the Renyi entropies are 
the most commonly used measures since they are sensi- 
tive to the long-distance quantum correlations of critical 
systems. 

In this paper, we study the Renyi entropy in the critical 
region of the anisotropic spin-s Heisenberg models for 
s = 1/2, 3/2, and 5/2. Consider an one-dimensional 
system of size L and composed by two subsystems A and 
B of sizes I and L — I, respectively. The Renyi entropy is 
defined as 



The first term, in this equation, is the conformal field 
theory (CFT) prediction in the scaling regime {L ^ I ^ 
1), and is given hy^''^^— 
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where pA is the reduced density matrix of the subsystem 
A. The von Neumann entropy is given by the limiting 
case a = 1. 

In the last few years a great effort has been made to 
understand the asymptotic behavior of Sa{L,l). It is 
expected that the ground state of critical one-dimensional 
systems gives a Renyi entropy that behaves as 



S^{L,l)^S'^''{L,l) + S°J-{L,l). 



for open boundary conditions (OBC). c 

charge, and cf^ and cf'^ are non-universal constants. 
The expression ([4]) for OBC has a small modification 

that is absorbed in the constant c{ \ as compared with 
the standard CFT expressioni^"— We chosed this ex- 
pression since, as shown in Ref. in the case of the 
XX model, it is accurate up to order (1/^) (see also Ref. 
[T3I ). It is interesting to mention that a generalization of 
the above equations for the excited states was proposed 
recentlyji^ 

The second term in Eq. Q seems to origin in the 
strong antiferromagnet correlations, as first argued by 
Laflorencie et al. in the case of the spin-1/2 Heisen- 
berg chain with OBCJ^ In open chains, since transla- 
tion invariance is broken, the energy density as a func- 
tion of the site I also shows a similar decaying alternat- 
ing term E°'"^. In Ref. [l^ by using bosonization tech- 
niques this oscillating term was calculated and compared 
with numerical evaluations of the entropy suggesting that 
(2) Sl'^'il) ~ E"'" - sin (7^)]^, where K is the Luttinger 
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parameter. In fact, these strong oscillations for the open 
chains were observed by several authors J^iHii^^^ 

More recently, Calabrese et al. in Ref. US (see also 
Ref. 113) investigate the anisotropic spin-1/2 Heisenberg 
chain with periodic boundary conditions and verified that 
even in this case those oscillations are still present if the 
Renyi index a >1. For the XX chain in a magnetic 
field with PBC/OBC, Calabrese and collaborators ob- 
tained exactly 5°**^, and observed an universal behavior 
of this oscillating termJ^i^i^i Based on these exact re- 
sults and on numerical calculations of the spin-1/2 XXZ 
chain with PBC at zero magnetic field, Calabrese and 
collaborators conjectured that S'^'^ has the following uni- 
versal behavior-i2iS^i21(see also Ref. [isl ) 
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we present some important relations that will be used to 
evaluate the Luttinger liquid parameter K. 

The ground state energy of a system of size L, as L ^■ 
oo, behaves as22i^ 
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where S — l{d ~ A) for the system with periodic (open) 
boundary condition, Vg is the sound velocity, Coo is the 
bulk ground state energy per site, and /oo is the surface 
free energy, that vanishes for the systems with PBC. 

The mass gap amplitudes of the finite-size corrections 
of the higher energy states, for a system with periodic 
(open) boundary conditions, are related to the anomalous 



dimensions 



'bulk 



(surface exponents x^). In the periodic 



case there are, for each primary operator Op {(3 — 1,2, ...) 
in the CFT, a tower of states ^ (L) in the spectrum of 
the Hamiltonian with asymptotic behavior— 



QOSC 



■ sm 



LP" 



{21 



1)4°) 



for OBC. The constants ga and g, 
and the critical exponents governing the decaying of the 
amplitudes of the oscillations are pa 
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are non-universal 
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In 



the PBC case [Eq. ([5])] the period of the oscillations 
depends on the Fermi momentum kp^ while in the OBC 



case [Eq. ^] it depends on kp = j^kp + 2{l+i) ■ 

The origin of the exponents pL^-* and , as observed 
by Cardy and Calabrese in Ref. [2^ are the conical 
spacial-time singularities produced in the conformal map- 
ping used to describe the reduced density matrix pA in 
the CFT. There are two important ingredients in the os- 
cillatory behavior of Eqs. ([5]) and the non- universal 
constants ga^ and g^"^ and the Luttinger parameter K 
(that gives the exponents pa'' and Pa^). Unfortunately, 
there is no prediction for g^'^ and gi") and as we are going 
to verify, in the case of the spin-s anisotropic Heisenberg 
model, they decrease dramatically as we increase the spin 
size s. 

Our aim in the present paper is to verify the above con- 
jectures more extensively for the spin-1/2 XXZ, as well 
for the critical regions of the spin s > 1/2 anisotropic 
Heisenberg model in the presence of an external z- 
magnetic field. The inclusion of the magnetic field is 
interesting because the magnetization and the Luttinger 
parameter K depend on its value. We have then more 
possibilities to verify the conjectures ([S]) and ©. The pa- 
rameter K in the case of the spin-1/2 XXZ chain can be 
calculated exactly from the Bethe ansatz solution of the 
model. For s > 1/2 the model is not exactly integrable 
but this parameter can be calculated numerically by ex- 
ploring the consequences of the conformal invariance of 
the quantum chain in the bulk limit. In the following 



eI^,{L) - Eo{L) = ^(xl,, + J + /) + o(L-i), (8) 

where ~ 0,1,2,.... For the chains with OBC the 
tower of states have energies^ 



E^{L) - Eo{L) = -^(.xf + j) + o{L-^l (9) 

with i = 0, 1,2, .... 

For models described by a Luttinger liquid CFT, which 
are the present cases, the Luttinger liquid parameter K 
is given hy K — where Xp is the lowest anomalous 
dimension obtained by using in Eq. ([5]) the lowest excited 
state in the sector whose magnetization is increased by 
one unit with respect to that of the ground state. 



II. THE MODEL 

We consider the anisotropic spin-s Heisenberg chain, 
also known as the spin-s XXZ chain, in the presence of a 
magnetic field h with Hamiltonian given by 



(10) 



were Sj — (s|, s^, s|) are the spin-s SU{2) operators and 
A = cos 7 is the anisotropy. 

We investigate the above model, using the density ma- 
trix renormalization group (DMRG)"^^ method with OBC 
and PBC, keeping up to m=4000 states per block in the 
final sweep. We have done ^6—11 sweeps, and the dis- 
carded weight was typically 10""^ - 10"^^ at that final 
sweep. In our DMRG procedure the center blocks are 
composed of (2s-|-l) states. 
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Let us first present some known results in the absence 
of the magnetic field, i.e, h = 0. It is well known that 
this model at the isotropic point A = 1 or 7 = is gap- 
less (gapful) for half-odd-integcr (integer) spins. '^^•■^•^ The 
anisotropic chains are critical and conformal invariant for 
— 1 < A < 1 with central charge c = 1 for half-odd- 
integer spinsj^"— On the other hand, in the case of in- 
teger spins a critical phase appears for —1 < A < Ac(s), 
where Ac(s) < 1 is a critical anisotropyi^ 

The spin- 1/2 XXZ chain is exactly soluble^i^ and for 
this reason some exact results are known on its critical 
region — 1 < A = cos 7 < 1. In particular, the anoma- 
lous dimension (surface exponent) associated to the low- 
est eigenenergy, in the sector with total spin ^-component 
= = 1' is given by Xp = ^ {x^ = 2xp) and 

the sound velocity is Vs = ZL™2^-42 j.-^^ ^^^^^ 

tion of the spin- 1/2 chain has also been explored in the 
context of the entanglement calculations!^"— Although 
exactly integrable, in this context, only some few issues 
were explored in the spin- 1/2 XXZ chain. This is due 
to the difficulty in extracting analytically results from its 
exact solution. 

For h ^ the model is in a critical and conformal 
invariant phase for he < \h\ < 1 + A. The critical field 
he = for |A| < 1, and for A >1, he = hc{A) depends 
continuously on A. In this phase the model is described 
by a Luttinger liquid phase whose parameter K — 

depends on the values of the magnetization per site m = 

m{h) = of the ground state, and the anisotropy A. 
The exact solution of the model allows us to obtain the 
exponents Xp — Xp{m). This is done by solving a set of 
non-linear integra l eq uations that, for the sake of brevity 
we refer to Ref. |4^. It is expected that for the half- 
odd spins s, with s >l/2, similar phases emerge when a 
magnetic field is applied. 



III. RESULTS 

As mentioned earlier, we need to calculate the anoma- 
lous dimension Xp in order to verify if the oscillating term 
of the Renyi entropies [Eqs. ([5]) and ([6])] decays with the 
new exponents pa^ = 2pa^ = For the spin- 1/2 case 
and h = 0, we know that Xp{m = 0) = For h ^ 

we can still determine Xp{m) exactly by solving numeri- 
cally a set of non-linear integral equations given in Ref. 
lisl . On the other hand, for s > 1/2 there are no exact 
results for Xp. 

In order to verify the dependence of the exponents pL^'' 
and with the Luttinger parameter K (or Xp), we need 
to use independent estimates for them. The independent 
estimates can be obtained from the mass gap of the eigen- 
spectrum and from the a— Renyi entropies. We are going 
to calculate these quantities using the DMRG technique. 
The estimate of Xp will be obtained by following basically 
the same procedure used by one of us in Ref. [l^. We 
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(0.25) 
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(0.2863) 
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0.3319 
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0.2963 




(0.575) 


(0.3334 ) 


( 0.3101) 


(0.2967) 


0.9980 


0.464 


0.358 


0.327 


0.309 




(0.49) 


(0.3599) 


(0.3264) 


(0.3086) 


2.0 




0.422 


0.363 








(0.4233) 


(0.3626) 





Table I: The anomalous dimension Xp for the spin-1/2 XXZ 
chain with PBC and some values of the magnetization m and 
anisotropy A. The values in between the parentheses are the 
exact ones (see text). 

consider the Hamiltonian defined in Eq. (|10|) with PBC 
to determine the anomalous dimension Xp. The value of 
Xp(m) is obtained from the limit L — > 00 of the finite-size 
sequences 

L[E{L,S^ + 1)-E{L,S^)] 
Xp(m,L) = — , (11) 

where as before = J2j ■ 

We estimated the sound velocity Vg using Eq. ([7]) with 
0=1. We also assume that Xp(m, L) behaves asymptot- 
icaly as 

Xp(m,L) = Xp(m) + ai/L"^ + a2/L'^, (12) 

where a; = — 4. These corrections are expected from 
the finite-size perturbation of the critical models (see Ref. 
[ssh . We use a simple fit procedure to obtain Xp(m) by 
considering typically system of sizes L ^ 16 — 96. 

A. s=l/2 

As a benchmark test, we consider first the spin-1/2 
XXZ chain. In Table I, we present the estimated values of 
Xp, for some values of magnetization m and anisotropy A, 
obtained from Eqs. (jlip and IT^ . We see on this table a 
clear agreement between our numerical results, obtained 
via DMRG, and the exact values shown in between the 
parentheses. We thus see that this procedure gives accu- 
rate estimates of Xp (similar results was found in Ref. [l^ 
for TO = 0). As mentioned before, for A > 1 the model 
is still critical in the region where the magnetic field pro- 
duces a non-zero magnetization (0 < to < 1/2).—"— Un- 
fortunately, for some values of A in this region, we ob- 
served that the DMRG is not stable. In particular, for 
systems with OBC and to > 3/10 we observed that the 
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DMRG is not stable for any A. For this reason, we were 
able to estimate Xp only for A < 2 and m < 3/10. 

Now, let us estimate the exponents pa^ and Pa\ Re- 
cently Calabrese et al. in Ref. [2j (see also Ref. |24I ) calcu- 
lated exactly these exponents for the spin- 1/2 XX chain 
(A = 0) with PBC. They also verified indirectly, through 
the DMRG algorithm, that pi^ = for A ^ and 

h = Q. Fagotti and Calabrese in Ref. also obtained 
p° exactly for the spin-1/2 XX chain in a magnetic field 
with OBC. We also include in our analysis the chains 
with OBC to show, indeed, that the conjecture holds for 
A 7^ 0. Here, instead of just verifying that the oscillating 
term 5'°'*'^ is consistent with a universal decay mediated 

by the exponents p^^'^ and as done in Ref. [23l , we 
are going to estimate directly these exponents. We obtain 
our estimates from the direct fit of the numerical data to 
the functions ^ and ([6]). 

In Figs. 1 and 2, we present the Renyi entropy Sa{L, I) 
as a function of I for the anisotropic spin-1/2 Heisenberg 
chains with PBC/OBC and for some values of m. The 
symbols (circle, squares, etc) are the numerical data and 
the solid lines [with the exception of Fig. 2(a)] connect 
the fitted points using Eq. ([2]) with c = 1. We show 
Sc,{L,l) only for I < L/2 since 3^(1, L- I) = Sa{L,l). 
As shown in Fig. 1(a), the amplitudes of the oscillations 
increase with the value of a. Similar results were also 
observed in Ref. [2^ It is important to stress that if 
a >1 we can only get a reasonable fit of the numerical 
data by considering the oscillating term 5'°'*'^, in addition 
to the standard term S^^'^ predicted by CFT. 

Note, from Eqs. ([5]) and that the periodicity of the 
oscillating term 5'°'*'^ depends on the value of the Fermi 
momentum kp. On the other hand, kp depends of the 
value of the magnetization, namely kp = (1/2 — m)7r. 
Therefore, if we change the magnetization, a change in 
the periodicity of the oscillations of Sa should be ob- 
served. According to Eqs. JS]) and dH), the period of 
oscillations is Al* = n/kp for PBC and Al* = n/kp 
for OBC. Our results presented in Figs. 1 and 2, are in 
perfect agreement with these predictions [see also Fig. 

md)]- 

In Fig. 2, we show the Renyi entropy for A > 1 
and some values of m. For |A| > 1 and m — the 
Renyi entropies tend to a constant j^i^ since the system 
is gapped. On the other hand, for non-zero magnetiza- 
tion (0 < m < 1/2) the system is critical, ^^"^^ therefore 
it is expected that Sa behaves as Eq. Indeed, we 

have observed these two behaviors, as illustrated in Figs. 
2(a)-(c). 

The numerical values of p^f ^ and pi""* obtained from 
the fitting together with the corresponding predicted ex- 
act values are presented in the Tables II and III. As 
we can see in these tables, the values found are very 
close to the expected ones (with the exception of A = 
cos(7r/50)=0.9980 and m = 0) which strongly suggest 
that the amplitudes of the oscillating terms decay as pre- 
dicted and observed previouslyii2ii^i2^i2i The fact that 
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Figure 1; (Color online). The Renyi entropy Sa{L,l) as a 
function of I for the spin-1/2 XXZ chain with A — \/2/2. (a) 
L = 96, PBC, m = and three values of a (see legend), (b) 
L = 180, OBC, m = 0, and two values of a (see legend). We 
added 0.3 in the values of the entropy for a = 2 in order to 
see both data in the same figure, (c) L = 96 with m = 1/6, 
and L = 80 with m — 3/10. In both bases a = 3 and PBC 
(d) L = 192 with m = 1/6, and L = 196 with m = 3/10. In 
both bases a = 3 and OBC. 
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Figure 2: (Color online). The Renyi entropy Sa{L,l) as a 
lunction of I for the spin- 1/2 XXZ chain for A > 1. (a) 
L = 96, A = 2, PBC, and three values of a (see legend), 
(b) S3{L,l) for chains with PBC, A = 2, and two values of 
magnetization m (see legend), (c) 5*3 (L,Z) for chains with 
OBC, A = 1.2 and two values of m (see legend). 



the estimated values of pa^ for A = 0.9980 and m = 
are not so accurate, is not a surprise since, at the isotropic 
point logarithmic corrections are present. Consequently 
the finite-size estimates have a very slow convergence for 
anisotropies close to A = 1. 



B. s>l/2 

Now, let us consider the case of spins s >l/2. As we 
will see below, it is not simple to confirm that the oscillat- 



ing term decays with the exponents pd 
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(0.3937) 




m — 


1/4 


0.675 

(0.6895) 
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Table II: The exponents p^2 ^ and p^^' obtained through a fit 
of Eq. m for the spin-1/2 XXZ chains with PBC, L = 96, 
and some values of A and m. We considered c = 1 in the 
fitting to Eq. ((2)). We also discarded the first points of Renyi 
entropy Sa in the fitting procedure (up to 5). The values in 
the parentheses are the exact ones. 
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to the small amplitude of the oscillations of S"^^, as we 



Table III: Same as Table II, but for chains with OBC. We fit 
the data with system sizes L = 180 — 196. We have discarded 
the first points (up to 12 points) in the fitting procedure. 



are going to see below. In order to get some idea of 
their order of the magnitude, let us consider the oscil- 
lations in the region around the middle of the chains. 
According to Eqs. ^ and (O the amplitudes are A^i'^ ~ 

5iP^(A,s)/LP^'"(^>-) and aI"^ ~ <?i°)(A, s)/LP^°'(^'^). 

(v) 

However, our numerical results indicate that gii and 
f;i°^are practically independent of the anisotropy value. 
The amplitudes A^J''^ and A^f^ are them expected to be 
smaller as p^a^ and pi""* decrease or, equivalently, as Xp 
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increases. For example, for A = \/2/2 (7 = 7r/4), the 
exponents are Xp = 0.333, 0.099, and 0.057 for s = 1/2, 
3/2 and 5/2, respectivelyi^^ Therefore, we expect, at this 
anisotropy, a decreasing of the amphtude of the oscil- 
lations as s increases. This can be observed when we 
compare Fig. [Ha) with the curves of Fig. ^a.) at 
A = Actually, our results indicate that for s — 3/2 

and s = 5/2 the amplitudes of the oscillations are ap- 
proximately one and two order of magnitude smaller, as 
compared with those of the s = 1/2 at A = respec- 
tively. With such small amplitudes it is quite difficult to 
extract, at the anisotropy A = \/2/2, accurate values of 
pi^-* and pa ^ for the spin chains with s > 1/2. 

In order to estimate the exponents p^f and pi"-* with 
some accuracy for quantum spin chains with s >l/2, we 
need to find at least a region in the critical phase of the 
model where Xp ~ 0.3, as in the spin-1/2 case. As it is 
well known, at the isotropic point A = 1 all spin-s Heisen- 
berg chains have Xp = 1/2. However, at this isotropic 
point, the operator governing the finite-size corrections is 
marginal producing logarithmic corrections that make a 
slow convergence in the finite-size estimates. We then se- 
lect an anisotropy close to the isotropic point. We choose 
A — cos7r/50 = 0.9980, although we should expect that 
even at this point the finite-size corrections are very large 
producing only rough estimates of Xp, p'a^ and 

In Table IV, we present the anomalous dimensions Xp 
for the spin-3/2 and 5/2 chains for A = 0.9980 and some 
values of magnetization m. Note that Xp <0.15, for the 
spin-5/2 quantum chains. This means that the evalua- 
tion of p'a'' and p'a' will be quite difficult even taking 
anisotropics close to the isotropic point. For this reason 
we concentrate in the s = 3/2 case at A = 0.9980 and 
m = 0, where the estimate value is Xp — 0.39. This value 
is close to the value Xp = 0.375 (see Table I) for the spin- 
1/2 with A = \/2/2. We then, naively, expect that the 
amplitudes A^if' and ^1°'* are of the same order in these 
two cases. However, to our surprise the amplitudes of 
the oscillations for s = 3/2 are quite smaller than those 
of the s = 1/2, as we can see in Figs. EJa) and[3{b). 
Similar results were also observed for the s = 3/2 chains 
with OBC and m 7^ (for a ~ 1 and m — 0, see also 
Ref. [22j). These results indicate that the non- universal 

constants g^ct' and gi""* appearing in Eqs. ([S]) and © are 
very small for s >l/2, as compared with the s = 1/2 

case. This means that the evaluation of pit' and p\t di- 
rectly from the decaying of the Renyi entropy oscillation, 
for the s >l/2 quantum chains, is a quite hard task. 

On despite of the above difficulties, we try to estimate 
pt and for the spin-3/2 chain at the anisotropy A = 
0.9980 and magnetization m — Q. As already mentioned, 
we should expect rough estimates since we are close the 
the isotropic point. Like the s = 1/2 case the oscillations, 
although smaller, happens only for a > 1. The central 
charge obtained by fitting Eq. ([2]) is, in the a = 1 case, 
c = 1.01 at A = ^/2/2 and c = 1.2 at A = 0.9980. If 



m = m = 1/6 m = 1/4 m = 3/10 

s = 3/2 0.39 0.1235 0.1177 0.1157 
s = 5/2 0.12 0.074 0.071 0.070 

Table IV: The anomalous dimension Xp for the anisotropic 
spin-3/2 and 5/2 Heisenberg chains with PBC, with 
anisotropy A — 0.9980 and some values of the magnetization 
m. 

we fix c = 1 the fitting is quite poor at A = 0.9980. A 
nice fit, at this anisotropy, is only obtained by allowing 
c as a free parameter. Certainly the estimate of c, as 
we increase L, decreases towards the exact value c = 1. 
Let us now concentrate on the case of a = 3 where the 
amplitudes are clearly present. We are going to fit the 
data with Eq. ([2]) in three distinct ways. 

Firstly, we try to fit the DMRG data to Eq. © by 
considering c = 1 fixed. In this case, the least square 
fitting give us g^f^ = for the best fit with = 0.0426^ 
[red triangles in Fig. ISKb)]. This means that through this 
procedure we get no oscillations^ [see Fig. ^h)] and we 
are not able to extract values of . 

In the second procedure we allow c and pg''' as free fit 
parameters. In this case we get ps = 1.21 and c = 1.27 
with = 0.0041 for L ^ 96 [black cicles in Fig. IlKb)]. 
Similar results were obtained for L = 48 where we get 

p(^) = 1.19.^ 

In the third procedure we fix the expected value of 
^ _I_ ^ 1/(6 X 0.39) = 0.427. In this case the get 

— 0.0087 which is higher than previous procedure. In 
Fig. EJb), the blue dashed line connects the fit to our data 

(p) 

by considering the predicted exponent pg — 0.427 fixed. 
Although the finite-size corrections are large, since we are 
close to the isotropic point, the results suggest that the 
oscillating term of the Renyi entropy decays as Eq. ([5|). 
We note that for the spin-3/2 chains, differently from 

the spin-1/2 case, the estimates of pa^ depend strongly 
on how many sites I of the subsystem A we discard. In 
FigEfc), we present the estimates of the exponents p'f ^ 

and pll'^ as function of l/l'^'"'"^, where Z'^**'^ is the number 
of sites we discard in the fit procedure (we discard the 
sites I = 1, . . . , Z*"'^). As we can note in this figure, for 

the spin-1/2 case the estimates of pa'' weakly depend 
on the values Z''**'^ and L. The weak dependence with 
the lattice size L also happens for s = 3/2. However, 
contrarily to the spin-1/2 case, those estimates are very 
sensitive to the number of sites l'^^'^'^ of the subsystem A 
we discard, as shown in Fig[3I^c). In order to take into 
account this effect, we assume that p'i^'\l'^'"=) behaves as 

p^HUn-pi^^ + i^^ + jj^. (13) 

If we fit our data for the spin-3/2 case [presented in 
Fig. |3Kc)] with this equation, we obtain p'f ^ = 0.49 and 
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Figure 3: (Color online), (a) The Renyi entropy Ss,{L,l) vs 
/ for the spin-3/2 and 5/2 chains of size L — 96, and two 
values of A (see legend), (b) 5*3 (L,Z) vs I for the spin-3/2 
chain of size L = 96 and A = 0.9980. The solid circles are the 
DMRG data and the other symbols are fits (see legend), (c) 
Estimates of p^*"' as function of l/l'^""^ for the spin-1/2 and 
3/2 chains with PBC. The symbols are the numerical data 
and the dashed lines connect the fitted data (see text), (d) 
Dio{L, I) vs I for the spin-3/2 chains at A = 0.9980 and two 
values of the magnetization (see legend). 



p'f^ = 0.34. These values are quite close to the predicted 
ones {p'f'' = 0.427 and pf^ = 0.32). In the above fit 
procedure the central charge c is also a free fit parame- 
ter. We also observed, in this case, that the estimates of 
c are sensitive to the values of and as we increase 
idisc ^YiQy get closer to the expected value c = 1. In 
FigMc), we present the fits for the spin-3/2 chains only 
for ^'^''"^ < 20, since for large values of I the amplitudes 
of the oscillations are of the same order of the numer- 
ical errors. Finding estimates of the exponent pa"^ for 
the spin-3/2 chains with OBC are even more difficult. In 
this case the estimates of are very sensitive to Z**^ and 
also to L. For this reason we are not able to find a sim- 
ple procedure to estimate this exponent for the spin-3/2 
chains with OBC. 

We could naively expect that for to > 0, where loga- 
rithmic corrections are not expected, pi^' would be better 
estimated. However, for m >0 the anomalous dimension 
Xp is small (see Table IV) and consequently the ampli- 
tudes of the oscillations are also small, complicating our 
analysis, as discussed earlier. However, we can observe 
an important feature of Eq. ([S]) in the Renyi entropy of 
the spin-3/2 chains. To better see this feature, it is con- 
venient to define the difference Da{L,l) — Sa — S^^^. 
In Fig. Efd), we present Diq{L,1) as a function of I for 
the spin-3/2 chain for A = 0.9980 and two values of 
the magnetization. We choose to present a large value 
of a since the amplitudes are bigger as we increase a. 
According to Eq. ([5]) the period of the oscillations is 
Al* = n/kp — 2/(2 — 2to). In fact, we have observed 
this periodicity, as shown in Fig. ^d) for two values of 
the magnetization. 



IV. DISCUSSION 

In this paper, we investigate the Renyi entropies of the 
spin-s anisotropic Heisenberg chains in a magnetic field. 
These quantum chains are critical and conformal invari- 
ant in a wide region of values of the magnetic field h 
and anisotropy A. The long-distance physics of this crit- 
ical region is described by a Luttinger Liquid CFT, with 
central charge c = 1. For this reason, these models are 
very attractive for testing predictions for one-dimensional 
critical systems. In particular, it is expected that the 
a— Renyi entropies have a term that oscillates with the 
subsystem size, whose amplitudes show a power-law de- 
cay with universal exponents p^^ and . These expo- 
nents are expected to depend on the Luttinger parame- 
ter K, i. e., pL^^ = 2K/a and = K/a for PBC and 
OBC, respectively. This universal behavior was obtained 
exactly for the spin-1/2 XXZ chains with a magnetic field 
with PBC and OBC for A = (XX chains) J^i^ More- 
over, for h = 0, DMRG calculations of the spin-1/2 XXZ 
chain also indicate that the oscillating term indeed de- 
cays as predicted. ^'^ As part of this work, we made an 
extensive study of the spin-1/2 model, but considered 
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a much wider region of couplings than those considered 
earher. Using DMRG technique, we also investigate ex- 
tensively the quantum spin-s chains (up to s — 5/2) with 
PBC/OBC for several values of anisotropy A and mag- 
netic field. Using the CFT machinery we were able to get 
accurate estimates of the Luttinger parameter K. For the 
spin-1/2 chains with PBC and OBC, we extract the ex- 
ponents pa^ and pa^ through a fit of Eq. ([2]) and confirm 
the predicted universal behavior of the Renyi entropy for 
several values of magnetization m and anisotropy A. For 
spin s = 3/2 our estimates of the the exponent pa are not 
so accurate due to the fact that the non-universal con- 
stants g^a^ and are very small (typically one order of 
magnitude smaller than the ones of spin-1/2 chain) . Even 

though, our results indicate that p^f' and pa ^ are related 
with the Luttinger parameter, as predicted. We also ob- 
serve that the periodicity of the oscillating term changes 
with the magnetization, as conjectured. For s > 3/2 

we were not able to extract the exponents pi^'' and pa ^ , 
mainly due to the fact that the ga ^ and ga ^ are so small, 
making the data difficult to analyze. In this case, the Lut- 
tinger parameter K is larger (compared with the ones of 



the spin-1/2) and also contributes to make the oscillat- 
ing term 5'°'"^ almost imperceptible. We could even think 
that this term is null as happens in the Ising model,— 
and that such small oscillations comes from numerical 
errors. However we are convinced that those small os- 
cillations are not related with the truncation errors in 
the DMRG. We have observed that those oscillations do 
not decrease as we increase the number states kept (up 
to m = 4000) in the DMRG procedure. These results 

(p) 

strongly indicate that the non-universal amplitudes ga 
and yi""* decrease very fast as the spin s increases. This 
makes a huge challenge the determination of the expo- 
nents pi^' and Pq"' with reasonable accuracy for s >l/2. 
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